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We study the dynamical melting of “hot” one-dimensional many-body localized systems. As dis¬ 
order is weakened below a critical value these non-thermal quantum glasses melt via a continuous 
dynamical phase transition into classical thermal liquids. By accounting for collective resonant tun¬ 
neling processes, we derive and numerically solve an effective model for such quantum-to-classical 
transitions and compute their universal critical properties. Notably, the classical thermal liquid 
exhibits a broad regime of anomalously slow sub-diffusive equilibration dynamics and energy trans¬ 
port. The subdiffusivc regime is characterized by a continuously evolving dynamical critical exponent 
that diverges with a universal power at the transition. Our approach elucidates the universal long¬ 
distance. low-energy scaling structure of many-body delocalization transitions in one dimension, in 
a way that is transparently connected to the underlying microscopic physics. 


I. INTRODUCTION 


The laws of thermodynamics can break down in dis¬ 
ordered quantum systems that are isolated from exter¬ 
nal heat sources, due to the localization of excitations 
that would ordinarily transport energy among distant re¬ 
gions to estabhsh thermal equilibrium [IHl]- Such many- 
body localized (MBL) systems have the remarkable prop¬ 
erty that all high-energy excited states behave like zero- 
temperature quantum ground states. Their spectrum 
can be labeled by an extensive set of local conserved 
quantities such that eigenstates exhibit boundary 

law scaling of entanglement entropy, characteristic of 
gapped quantum ground states. This raises the intrigu¬ 
ing possibility that quantum coherent phenomena, typi¬ 
cally associated with zero-teraperatme systems, can oc¬ 
cur in arbitrarily “hot” matter. Examples of such quan¬ 
tum coherent phenomena include symmetry breaking be¬ 
low the equilibrium lower critical dunension [5], topo¬ 
logical edge states [51ll2|. and quantum criticality |13h 
I15| . MBL systems also host novel out-of-equilibrium dy¬ 
namical phase transitions, where thermodynamics breaks 
down sharply at a critical point |1()H19J . Cold atomic, 
molecular, and trapped ion systems offer a promising ex¬ 
perimental platform to explore these theoretical ideas. 
Indeed, issues of thermalization and excited state dynam¬ 
ics necessarily aidse in such systems, as they are inher¬ 
ently well isolated from their surroundings, and typically 
camiot be cooled to low temperatures (compared to their 
characteristic energy scales). 

Theoretical investigations uito these questions, how¬ 
ever. must tackle a daunting combination of out-of- 
equilibrium quantum dynamics, interactions, and disor¬ 
der. Consequently, most existing theoretical work on 
MBL systems has been confined to small-scale numer¬ 
ics and analysis of phenomenological models. Remark¬ 
ably, because of the short localization length and short- 
range entanglement structure of MBL eigenstates deep 
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FIG. 1. Phase diagram and finite-size crossovers for 
onc-dimcnsional many-body delocalization transitions. W pa¬ 
rameterizes disorder strength. For strong disorder W > Wc. 
a non-crgodic many-body localized glass is obtained, which 
lacks thermal transport but exhibits slow dephasing and en¬ 
tanglement dynamics with length-time scaling L ~ log Tdsph 
(where the subscript “deph” emphasizes that this is the dc- 
phasing time associated with virtual quantum fluctuations 
rather than the energy transport time scale, which is strictly 
infinite in the strong disorder glass). Below a critical disorder 
strength. Wc- the MBL glass melts into a thermal, energy¬ 
conducting liquid. The melting transition is continuous (sec¬ 
ond order), characterized by a single diverging length-scale 
^ ~ \W — Wc\~'', where i/ w 3.5 ± 0.3. Energy transport (en¬ 
tanglement spreading) in the thermal liquid arc characterized 
by a subdiffusive (sub-ballistic) power law scaling between 
time and length r ~ with dynamical exponent 2 ~ ^ that 
diverges continuously upon approaching the transition. 


within the locahzed phase, these approache,s have met 
with considerable success in gleaning properties of the 
MBL phase itself. In c:oiitrast, these same techniques are 
poorly suited to access the universal properties of disor¬ 
dered criticality and dynamical phase transitions out of 
the locahzed phase. 

For example, whereas the ground states of disordered 
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oiie-dimensional systems are always localized, the ex¬ 
cited states are only many-body localized for sufficiently 
strong disorder, and melt into a self-thermalizing in¬ 
coherent classical fluid at weaker disorder. These two 
dynamical phases are separated by an apparently con¬ 
tinuous many-body delocalization phase transition [3]. 
Such excited-state delocalization transitions are neither 
classical thermal phase transitions nor zero-temperature 
quantum phase transitions, but rather, represent a novel 
class of dynamical qnantum-to-classical criticality. Such 
transitions are driven by long-distance properties whose 
characteristic length scale diverges as the critical point is 
approached, rendering traditional numerical methods in¬ 
creasingly nmeliable near the critical point. For example, 
recent numerical studies |lf^) of a many-body delocaliza¬ 
tion transition obtain critical exponents that contradict 
fundamental bounds j2()) . suggesting that small scale nu¬ 
merics cannot access the true long-distance scaling prop¬ 
erties. New approaches are clearly needed to understand 
the nature of this new class of dynamical transitions. 

To this end, we develop an effective model for the for¬ 
mation of collective many-body resonances that destabi¬ 
lize the MBL phase at weak disorder, in order to compute 
the universal scaling properties of many-body delocaliza¬ 
tion transitions. While the problem of identifying generic 
many-body resonances is as hard as directly solving the 
full quantum problem, at any continuous phase transi¬ 
tion, one expects a self-similar hierarc'hical scaling struc- 
tme. This suggests that the critical resonant cluster may 
be constructed hierarchically in the spirit of the renor¬ 
malization group. We implement an efficient numerical 
procedure to identify the formation of such delocalizing 
resonant clusters and study the critical properties of the 
resulting effective model. Our construction is motivated 
in part by the requirement that the MBL transition it¬ 
self must be thermal |2l) and thus, incoherent and clas¬ 
sical. Accordingly, the critical properties should be well 
describe<l by an effective classical statistical mechanics 
model, e.g. which ignores quantum interference effects 
such as those responsible for weak-localization of non¬ 
interacting fermions in equilibrium settings. 

We obtam scaling results for the critical properties 
of the disorder-strength-tuned MBL transition in one- 
dimension. by numerical simulation of the resonance 
model. We find a continuous dynamical phase transition 
characterized by a diverging length scale fy in which a 
quantum MBL glass melts into a classical thermal liquid. 
The predicted phase diagram and cross-over structure of 
this effective model are summarized in Fig. 

The effective resonance percolation model also pro¬ 
vides insight into the dynamics of thermal transport and 
entanglement of the deloc'alized thermal liqui<l, near the 
melting transition (see Fig.[^. In one dimension, we find 
that transport on the deloc’alized side of the transition is 
anomalously slow, characterized by power-law subdiffu¬ 
sion with a continuously evolving dynamical exponent, 2 , 
which diverges in a universal fashion at the MBL tran¬ 
sition. We argue that this subdiffusion stems from tun¬ 


neling through rare insulating regions, in agreement with 
recent numerical studies jl6lll7| and a phenomenological 
renormalization group approach |19) . and establish a scal¬ 
ing relation between the correlation length exponent j/, 
characterizing the divergence of the many-body localiza¬ 
tion length fy and the exponent governing the divergence 
of 2 . 

Tills ID MBL transition was recently studied by mean 
field methods [22], exact diagonalization [IMT51I22] and a 
phenomenological renormalization group procedure |19j . 
Exact diagonalization studies on 10 — 20 site chains ob¬ 
serve an apparently continuous transition, but obtain a 
correlation length exponent that is inc'ompatible with 
fundamental bounds j20| . suggesting that the observed 
scaling properties are likely finite size artifacts that do 
not reflect the true long-distance universal properties of 
the MBL transition. In contrast, the critical exponents 
obtained in our approac'h comply with the Harris-Chayes 
bounds [20] . 

Onr results are broadly consistent with those of 
Ref. 1191 where available. We also comment on the many- 
body delocalization transition of melting quantum criti¬ 
cal glasses (QCGs) [T5HT^ [21) . critical analogs of MBL 
phases that arise, for example, in random-bond spin 
chains. We argue that the delocalization transition for 
these QCG phases has identical properties to the ordi¬ 
nary ID MBL transition. 


II. MODEL AND RESONANCES 


Though onr effective model for the MBL transition is 
largely insensitive to microscopic choice of a particular 
model, for concreteness, we will phrase our discussion in 
terms of a paradigmatic MBL model, the random-field 
XXZ chain: 

77mbl - . 

( 1 ) 


Here, are Pauli-matrices for spm-I/2 degrees of 

freedom. This XXZ chain maps, via a standard .lordan- 
Wigner transformation, into a chain of spinless interact¬ 
ing fermions, where .7 is the hopping amplitude, V J 
is the interaction strength, and /ij is the chemical poten¬ 
tial uniformly distributed on the interval [0,W], where 
W serves as a measme of disorder, and we use the spin 
and fermionic languages interchangeably. We will work 
in the Fock basis of single-particle orbitals, (pi, of the 
non-interacting part (y = 0) of Q, which have en¬ 
ergy ei « /it -|- C7 ^re localized near site i. 

We will then consider the interactions 17 as a perturba¬ 
tion expressed in this basis, though we remark that our 
method is not confined to low-order perturbation theory, 
but rather includes contributions from all orders in V. 
We will be interested in the properties of typical highly 
excited states in the middle of the many-body spectrum. 
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This Haiihltoiiian has been considered as a prototypical 
example of a many-body localized (MBL) system (see 
e.g. |25| for a recent review), characterized by the ab¬ 
sence of thermal transport and by a slow, logarithmic 
spreading of entanglement due to dephasing 0 1261127) . 
When F / 0. the system is beheved to exhibit a MBL- 
to-thermal dynamical phase transition as a function of 
disorder strength |l6H19| . which separates a iion-ergodic 
MBL quantum glass phase from a ergodic high tempera¬ 
ture classical liquid. 

To motivate the problem of delocalization, imagine cre¬ 
ating an energy wavepacket initially localized near a sin¬ 
gle site in the MBL system. How does such a packet 
evolve with time? Deep in the localized phase, it is 
extremely unlikely that it can tuimel from a given site 
(bond) to another, due to the wide variance of site (bond) 
energies. More precisely, the amplitude to tunnel through 
a distance x of the localized phase is 

r(rr) « (2) 

where itq In , is the single parti¬ 

cle localization length (measured in units of the lattice 
spacing), where we have chosen an expression that in¬ 
terpolates smoothly between the strong disorder limit 
xo ^ hi”^ and weak disorder limit xq ~ {^) • 
If the tunneling strength, r(xij), between spin pairs i 
and j separated by distance Xjj is much smaller than 
their energy difference SEij, then the true many-body 
eigenstate is very close to a product state of indepen¬ 
dent configurations of each bond. On the other hand, if 
r(xij) > SEij, then the many-body eigenstate contains 
an entangled superposition of the degrees of freedom on 
the two sites, and we say that the two sites are resonantly 
linked. We estimate the probability that a given site is 
resonantly linked to at least one other is (1 — e”"^”), with 
Ar = zy(0) dxr{x) « ^3;o, where i/(0) » 1/J is the 
single particle density of states. 

For very strong disorder (xq <gC 1), the density of reso¬ 
nantly linked pairs is very small, prcs ~ Ar? 1, so that 
the resonant links are well isolated from each other, and 
they do not disrupt the properties of the surrounding lo¬ 
calized phase. As disorder is weakened, these resonant 
links occur more frequently, and eventually disrupt lo¬ 
calization below a critical disorder strength Wc Since, 
in a thermal phase, all spins are highly inter-entangled, 
delocalization and the recovery of thermalization at weak 
disorder must occur via a highly collective resonance in¬ 
volving every single spin. Identifying ail possible generic 
many-body resonances is a hard problem, and in fact as 
difficult as exactly solving the full quantum many-body 
problem. However, near the delocalization critical point, 
it is natural to expect that the critical resonant cluster ex¬ 
hibits scaling structme that is geometrically self-similar 
(in a statistical sense) at different energy scales. An¬ 
ticipating this hierarchical structure motivates a simpler 
way of identifying collective resonances by first identify¬ 
ing small resonant clusters (e.g. resonant pairs as above). 


and then examining whether groups of these small res¬ 
onant clusters can collectively resonate on some lower 
energy scale (longer timescale), and so on. Using this 
perspective, we will construed an effective quantmn per¬ 
colation model for eielocalization by hierarchical collec¬ 
tive resonances. 


III. EFFECTIVE MODEL FOR THE 
DELOCALIZATION TRANSITION 

A cluster of m resonantly linked sites is characterized 
by bands of A*» 2"’ many-body energy levels, and thus 
has an exponentially small level spacing in m. Roughly 
speaking, collective reai'rangements among spins of two 
such clusters containing mi ^2 sites respectively and sep¬ 
arated by distance L 12 , can occur resonantly if the effec¬ 
tive tunneling T 12 is larger than the level spacing of the 
combined cluster, scaling as ~ We therefore 

see that this exponential suppression of the level spacing 
will favor large, many-spin resonant clusters. This moti¬ 
vates an effective model for the delocalization transition, 
consisting of the following iterative procedure: 

1. Examine a disordered chain to identify resonantly 
linked pairs of spins, and group them into 2-spin 
clusters 

2. Identify the strongest resonantly coupled pair of 
clusters (i.e. the one with largest interaction Fij) 
and merge them into a new larger resonant cluster 

3. Account for the reduction of level spacing on newly 
formed resonant cluster 

4. Compute the effective interaction of the new cluster 
with other clusters (as described in detail below) 

5. Examine whether the reduced level spacing enables 
any of the newly formed clusters to resonantly in¬ 
teract, if so merge these clusters 

Steps 2-5 can then be iterated (numerically) until all pos¬ 
sible clusters are formed. 

Before giving a detailed account of this procedure, 
we briefly sketch its potential outcomes and underlying 
philosophy. For strong disorder, the resonance merging 
ceases after a small number of steps, producing only di¬ 
lute well isolated small resonant clusters embedded inside 
strongly localized quantum glass. Upon weakening dis¬ 
order, the size of the resonant clusters increases, diverg¬ 
ing continuously at a critical value of disorder W = W^, 
where a single cluster just barely percolates across the 
system (in the limit of infinite system size). For weaker 
disorder (H^ < Wc), the resonance merging proc:edure 
continues until all degrees of freedom are subsumed into 
a single percolating cluster (in the infinite system limit), 
resulting in a delocalized thermal phase. Nevertheless, 
we will see that the near critical transport properties of 
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the delocalized phase exhibit anomalously slow dynam¬ 
ics, with an unusual scaling structure, due to the prox¬ 
imity to the quantum glass phase. 

In this procedure, resonances arise from small mis¬ 
match in energy scales, corresponding to a slow beating 
frequency or long time in dynamical processes. Over the 
coarse of iterations, the dynamical time scales associated 
with the merging of resonances steadily lengthen. In this 
way, the effective model closely resembles a renormaliza¬ 
tion group procedure - first accounting for fast degrees 
of freedom and progressively proceeding to slower and 
slower modes. We also remark that, in principle, the 
problem of identifying generic many-spin resonances is 
factorially complicated in system size, and hence is no 
simpler than directly solving the full quantum problem. 
From our analytic treatment of the strong disorder phase, 
we expect that the most likely type of resonances are 
those between two bonds (orbitals), motivating our use 
of step 1) above to “seed” the iterative resonance merging 
procedure. Moreover, near the transition, whether or not 
a resonant cluster percolates across the system is driven 
by highly collective resonances among large clusters of 
many spins, and hence the low-energy, long time scale 
dynamics of the near-critical cluster are expected to be 
insensitive to the details of early steps in the procedure 
(i.e. are universal). 


A. Criterion for merging resonsint clusters 


To test if two clusters i and j can resonate (step 2), we 
compare their coupling Fij to the energy mismatch, SEi^j 
of their energy levels. For two large clusters i.j, such 
that the bandwidth of each cluster exceeds the level spac¬ 
ing of the other (min{Ai,Aj} > max{(5i, (fj}), we define 
an energy mismatch between two clusters i,j as: 

■ f*. , 1 . In the other case, where one cluster's band- 


width is smaller than the other’s level spacing, we define 
5Ei j = max{^i,(5j} — min{Ai, A^} (see Appendix A1 1 . 
We take the bandwidth of the newly merged cluster to 
be: Aiuj = Aj -f Aj + Fjj, and the level spacing to be 


6i\jj = —-■ We emphasize that the main ingredi¬ 

ent at this step is the exj>onential reduction of the level 
spacing when two clusters are merged: we ciiecked nu¬ 
merically that the universal properties of the transition 
are insensitive to the other details of the procedure. 


B. R.enormalization of inter-cluster couplings 


After merging two clusters, one should in principle 
compute a new set of inter-cluster couplings (step 4 
above). Consider merging two clusters i and j, with 
number of spins mi-rrij and bandwidths Ai,Aj respec¬ 
tively, and interacting with coupling Fj^j. In order for 
a third cluster, k. to drive a c:ollective rearrangement of 
all levels of the merged cluster, i Uj, and hence 


access the full many-body level spacing of i U k must 
interact separately with both i and j. The effective inter¬ 
action strength for such a process depends on the energy 
mismatch between k and i,j respectively, and can oc¬ 
cur by several different possible processes. For instance, 
consider the process k i j, where k first interacts 
with i and then i with j. If the energy mismatch, SEk,i/j 
between k and i or j is large compared to the couplings 
between k and i or j , then the interaction between 
k and iU j can be computed pertiirbatively: 

~ 6E,,i ■ 

In the alternative case where all three clusters i, j, and 
k can resonate , e.g. F^^j > 5Ek,i, we may estimate 
the effective c:oupling' between k and iU j by classically 
adding the time taken for an excitation to transfer from 
fc —> i —> J: 

+ (4) 

There are a few possible routes for k to excite both i and 
j - for example k can first excite i, then i excites j, or 
k can excite both i and j, etc. We compute F,t,iuj for 
each process using the ‘perturbative’ or ‘classical’ rules 
above as appropriate, and c:hoose their maximum as the 
new effective couphng between k and i U j. 

Computing the renormalized inter-cluster couplings in 
step 4 is conceptually important to avoid a potential in¬ 
stability of the strong disorder phase (see Appendix [Aji. 
however, in practice we find that this instability is absent 
for all numerically accessible system sizes and find identi¬ 
cal scaling independent of whether or not the inter-chister 
couplings are renormalized. We remark that the way 
the couplings are renormalized is remuiiscent of the ap¬ 
proach of Ref. jl9) where metallic and isolating bloc:ks are 
merged by using a renormalization group scheme. How¬ 
ever, we point out that we only merge “metallic” (reso¬ 
nant) clusters in our model, thus avoiding the thorny is¬ 
sue of determining the proper renormalized couplings af¬ 
ter merging “mixed” insulating and thernializmg blocks. 


IV. SCALING STRUCTURE OF THE MBL 
TRANSITION 

To identify the scaling structure of the transition pre¬ 
dicted by the effective resonance percolation model de¬ 
scribed above, we examine the dependence of the size, 
^loci (number of spins) of the largest resonant cluster as 
a fnnc:tion of disorder strength, W, for various system 
sizes, L. Disorder averaged results for V = 0.3 are shown 
in Fig. For W ^ Wc one obtains a percolating cluster 
with ^loc ~ L with probability close to one. In the fi¬ 
nite size systems simulated, ^loc/L crosses over from one 
(If ' <gC Wc) to zero (W Wc) as disorder is increased. 

The curves for different system sizes cross at the same 
value of disorder strength. whic:h we identify as the crit- 
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FIG. 2. Finite-size scaling for the disorder averaged 
localization length ^loc (defined as the size of the longest 
resonant cluster for a given disorder realization). col¬ 

lapse to a universal form (main panel) when plotted against 
the scaling variable {\V — , with correlation length 

exponent v w 3.5 ± 0.3. Insets show the unsealed curves. 
Results arc averaged over ~ 10"* disorder realizations. 


ical disorder strength Wc- Moreover, the data for dif¬ 
ferent system sizes collapse to a universal scaling form 
^loc = L =. {(Td^ — with correlation length ex¬ 

ponent V = 3.5 ih 0.3, as shown by plotting against 
the scaling variable {W — Wc)L}^^ (Fig- Here E(a;) 
is a universal scaling function interpolating between one 
(x <C 0) and zero (x 0). 

This scaling collapse indicates the existence of a length 
scale ^ that diverges when the transition is approached 
from either side as ^ ~ \W — 1141“", signaling a sharp 
continuous delocalization phase transition. The diverg¬ 
ing length scale ^ can be interpreted as the localization 
length for W > Wc. and as we will see below, character¬ 
izes the length of insulating gaps in the transport path 
for W < Wc. 

We remark that the value of the scaling fmiction at 
criticality: E(0) 10“^, is anomalously low compared 

to ordinary percolation (for which H(0) would be fts 0.5). 
This strong asymmetry indicates that the transition is 
driven by rare resonant clusters that are quite sparse and 
widely separated (e.g. clusters of size i are separated by 
typical distance E(0)“^ x ^ ^ at criticality). 

We also performed simulations of a simpler model that 
ignores the renormalization of the couplings when merg¬ 
ing clusters. As claimed above, despite the potential sta¬ 
bility issues of the strong disorder phase within that sim¬ 
plified model, we find essentially identical universal per¬ 
colation curves using these simplified rules, compatible 
with the same imiversal exponent u « 3.5 within error 
bars (see Appendix [A)|. 


V. NEAR-CRITICAL DYNAMICS IN THE 
DELOCALIZED PHASE 


We now turn to the task of computing the near-critical 
scaling of transport and entanglement dynamics, in the 
delocalized phase (W < Wc). Our goal will be to mi- 
derstand the scaling structure underpinning the critical 
slowing down of energy transport as W approaches Wc 
from below. Since the delocalized phase is thermal and 
at high temperature |21| . transport occurs via the ther¬ 
mally incoherent transfer of energy among bonds in the 
resonant cluster. Hence, we model the spread of excess 
energy (initially localized near position x as a function 
of time) as a classical random walk across the resonant 
cluster, with a timescale tab — ^1 J{Lab) to transfer 
excitations between resonantly coupled bonds A and B 
separated by distance Lab- 

The iterative merging procedure yields detailed infor¬ 
mation about the connectivity structure and the coupling 
strengths for each participating link in the resonant clns- 
ter. Here, it is important to note that, while all degrees 
of freedom are part of the percolating resonant cluster 
for W < Wc, not all spins contribute equally to trans¬ 
port and dynamics. Rather, dynamical properties are 
dominated by a subset of efficiently connected spins that 
form the “backbone” of the transport path thi'ough the 
resonant cluster. Then, using standard Green’s function 
methods, we compute the time evolution of an initially 
well-localized energy wave packet spreading across the 
resonances network via a random walk. Our main re¬ 
sult is that the delocalized phase exhibits a broad regime 
of anomalously slow sub-diffusive equilibration dynam¬ 
ics and energy transport. We observe tliis subdiffusive 
behavior in our model using both renormalized ami non- 
renormalized couplings, but find that finite-size effects 
are mucli weaker for the simplified model that ignores 
renormalization, thus allowing us to extract much cleaner 
data. We will therefore focus on this simplified model in 
the following - recall that it reproduces equally well the 
critical percolation curves described in the previous sec¬ 
tion p5] . 

Averaging over initial position x, and over disorder 
configiuations, we find that the mean-square displace¬ 
ment of an excitation grows as a power law in time 

{Fig.[^, 

(5) 


For ordinary classical diffiLsion, z = 2. In contrast, we 
find that z{W) diverges in the glassy phase, indicating 
an absence of energy transport and a breakdown of ther¬ 
mal equilibration. On the delocalized side of the critical 
point, l/z(W) increases continuously from zero as a func¬ 
tion of detuning from the critical disorder strength: 


z(W) 


ZQ 

{Wc - wf ’ 


(W< W,). 


( 6 ) 


Though seemingly a distinct universal exponent, Q is re¬ 
lated by a general scaling relation to the log-dynamical 
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FIG. 3. Energy transport near the trEinsition. In the 
delocalized phase {\V < W^) excess energy initially localized 
near the origin spreads to distance St ~ in time t (bottom 
inset). The dynamical critical exponent z diverges continu- 
oush' from the delocalized side of the transition, and vanishes 
inside the quantum critical glass phase (main panel). The nu¬ 
merical results are consistent with z diverging as \W — \Vc\~‘’ 
w'ith (^ = i' (top inset). 


exponent ij} of the localized phase, and the correlation 
length exponent u. To see this, note that at the critical 
point, transport along the critical resonant chain occurs 
by tnniieliiig through arbitrarily long localized regions, 
hence energy scales with length as ~ e~^ at the critical 
point, just as in the localized phase. This must cross over 
smoothly to the power law scaling of on length scales 
L « implying the scaling relation: 

C - (7) 

Using this scaling relation and the dehnition of the cor¬ 
relation length exponent we may rewrite Q as z ~ ^ 
(ignoring logarithmic corrections and non-universal pref¬ 
actors). Accessing the scaling regime Q to measure 
numerically is difficidt because of the form of the finite- 
size percolation emves (Fig.[^: close to the critical point 
with W < Wc, ^loc > L even for moderately large sys¬ 
tems. In particular, the small non-zero value of z~^ at 
Wc is a finite size artifact, and upon subtracting the value 
of z[Wc)~^, we find critical behavior compatible with the 
expected scaling relation C, = v (Fig. [^upper inset). 


A. Structure of transport path 

Subdiffusive transport arises from the broad distri¬ 
bution of effective tunneling finks, Fij, in the resonant 
backbone, corresponding to a power-law distribution of 


timescales |171 12[j| . Tij ~ l/Fiji 

with 1 < Q < 2 such that the mean waiting time 
/o '''p{T)dT is divergent. To see the relation between this 
power-law distribution and subdiffusive transport prop¬ 
erties, note that energy transport occurs as a random 
walk on a one dimensional chain whose '‘sites’’ are two- 
spin bonds in the resonant cluster, and whose links are 
weighted by a waiting time r corresponding to the in¬ 
verse effective coupling between “sites” (In principle the 
bonds also have variable lengths, ^ ~ log r. However, 
unlike the waiting-times the bond-lengths have a finite 
mean, i, ami in the following we approximate the length 
of each step as £). In N steps, a random walker traverses 
~ 0{\/N) different bonds, revisiting each approximately 
yfN times, and moving total distance L w '/N (in units 
of the average step size). For broad distribution of wait¬ 
ing times, the time taken for N steps is dominated by 
the longest waiting time, T, encomitered. The probabil¬ 
ity of encoimtering a bond with waiting time of at least 
T is P{t > T) = p(t) ss We can reasonably 

expect to find a link with waiting time t > T among 
L =3 y/N different bonds if LP{t > T) « 1, and there¬ 
fore the slowest link in region of size L has waiting time 
T{L) S3 . Smee this weak link is revisited order 

\/N S3 L times, the total time to move distance L scales 
as tonergy(7>) ^ . Comparing to (j^, we identify the 

dynamical exponent for energy transport as: 

Q 

■2 = ^energy = Ti (9) 

a — 1 

implying that a approaches one near the transition as: 
q(W)-1~|W-Wc|^. 

Note also that the broad distribution of waiting 
times (j^ corresponds to a probability of encountering 
a link of length £ = log r in the transport path along the 
resonant backbone, near the transition, given by 

= ( 10 ) 

Hence the typical spacing between “vertebrae” in the res¬ 
onant backbone is ss corresponding to energy scale e“^. 
However, in a region of length L 3> C it is extremely 
unlikely to avoid encountering an atypically long gap of 
length f*, defined by jP{£*) 1, i.e. ^ ^log 

Such long rare links involve waiting time r ~ e^' ~ L^, 
which dominates the time required to traverse a segment 
of size L. 

B. Entanglement dynzunics vs. thermEil transport 

We note that the dynamics of energy transport scale 
differently than those of entanglement since energy is a 
conserved quantity that can oidy be transferred among 
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subregions whereas entanglement can be freely gener¬ 
ated |191 130} ■ Hence, whereas energy transport can be 
viewed as a random walk of conserved particle-like excita¬ 
tions. entropy spreads deterministically in all directions 
simultaneously j30) . In order to entangle regions sepa¬ 
rated by distance L , entanglement must spread across 
order L links, only visiting each once (in the same spirit 
as the second law of thermodynamics, once a bond is 
entangled with many others, it is extremely unlikely to 
later disentangle itself). Therefore the timescale for en¬ 
tanglement spread is dominated by the longest typical 
waiting time encountered, tent ^ T(L) fts L"-‘. Com¬ 
paring, we find a different effective dynamical exponent 
for the spread of entanglement: 


2ent 



-energy 


- 1 , 


( 11 ) 


with ~ 2 ~ C Hear the transition. 


C. Seeding of optical conductivity 

So far we discussed transport from the perspective 
of the propagation of an initially well-localized energy 
wave packet. This is natural in ultracold atomic sys¬ 
tems (see e.g. j31)l. The power-law distribution of time 
scales also implies an anomalous power-law frequency 
dependence for AC conductivities. For systems with a 
conserved spin-component or particle number, one can 
compute the scaling of the optical number conductivity 
via the Einstem relation ~ Z)(w)y, where D{aj) is 
the (frequency-scale dependent) diffusion constant, and 
X is the static compressibility (which will be constant 
throughout the phase diagram). The diffusion c:onstant 
scales like D ~ uL'^ ~ where z = Zenergy is the dy¬ 

namical scaling of energy, hence the optical conductivity 
will also scale as: 

( 12 ) 

This expression interpolates between constant conduc¬ 
tivity in the diffusive limit (z = 2) for weak disorder, 
to cr{Wc) ~ w at the MBL transition z = oo. This scal¬ 
ing argument agrees with the recent rare-regions analysis 
and mmierical results of Ref. HZ]. 

In two dimensions, we expect a non-zero diffusion con¬ 
stant and Bnite DC conductivity in the entire region 
W < Wc- By the above argument, our effective resonance 
model again predicts z = co (tunneling-like) dynamics 
at the MBL transition in higher dimensions, suggesting 
ct(cj) ~ uj in the critical regime. 


VI. SUMMARY OF RESULTS 
A. Phase diagram 

Let us smnmarize our results in the thermal phase 
(IT < Wc)- Near the critical point, the system exhibits 


anomalously slow subdiffusive dynamics t L^. Inside 
the critical regime (i.e. on length scales L \W — 
lTc|~‘'), the critical resonant cliLSter mediates coherent 
transport of energy with characteristic scaling of length 
and time logr ~ L. This scahng relation is reminisc:ent 
of the one governing the dynamics of dephasing and en¬ 
tanglement growth in the localized MBL phase, with the 
important distinction that dynamics in the strong disor¬ 
der MBL phase describe only virtual (“’off-shell”) transi¬ 
tions rather than real (“‘on-shell”, or resonant) processes 
required for thermal transport. Heiic:e, whereas the MBL 
glass is non-ergodic:, the energy eigenstates at the glass 
melting critic:al point are ergodic and thermal, with vol¬ 
ume law scaling of entanglement in every eigenstate (con¬ 
sistent with general entanglement monotonicity reqnire- 
ments [21] 1. However, while the pure eigenstates at the 
delocalization transition exliibit thermal behavior, start¬ 
ing from a superposition of energy eigenstates (the only 
initial conditions that can be prepared experimentally in 
finite time), the system will take super-polynomially long 
time (in system size) to equihbrate and thermalize. The 
resulting phase diagram and crossover scales are shown 
in Fig. 


B. Quantum critical glasses and extended 
interactions 


We now briefly comment on the generalization of our 
results to critical analogs of MBL systems, dubbed quan¬ 
tum critical glasses (QCGs) in Ref. 1131 Examples of 
QCGs include phase transitions between MBL phases 
with different (discrete) symmetry or topological or¬ 
der [l5IfT^ [211, as well as stable excited state critical 
phases |15j . At strong disorder, QCGs are characterized 
by a stretched exponential tunneling in which energy, E, 
scales with distance, a:, as: £ ~ with universal 

exponent ip < 1 whose value depends on the universal¬ 
ity class of the QCG in question. As for the MBL case, 
upon weakening disorder, one expects a continuous de¬ 
localization transition from QCG to thermal liquid in 
which resonant clnsters percolate. 

Whereas the renormalization of inter-cluster couplings 
plays little role in the scaling properties of the MBL tran¬ 
sition, we expect it may play a more crucial role in the 
QCG delocalization transition. Loosely, two clusters with 
number of spuis n \^2 respectively, can inter-resoiiate if 
their coupling ri 2 exceeds the many-body level spacing 
on the clusters, w As the clusters are formed, 

the bare coupling between them, ro(xi 2 ) ~ 

(where x \2 is the inter-clnster distanc:e), is perturba- 

tively reduced by a factor of ~ 

where oj^yp is the typical spacing of spins in the clus¬ 
ters, and Atyp is their typical energy mismatch. At 
the delocahzation transition, one expects the clnsters 
to just barely resonate, i.e. ni -|- 7i2 *3 x\ 2 > giving 






effective inter-cluster coupling ri 2 « e ■^^^■^^ro(3;i2), 
where « (xtyp/xo)'^ — logAtyp. As the exponen¬ 
tial prefactor dominates over the stretched exponential 
bare coupling To, this indicates the renormalization of 
couplings gives effective exponential-in-distance inter¬ 
actions between large resonant clusters in QCGs near 
their delocalization transition, in contrast to their longer 
range stretched-exponential-in-distance interactions be¬ 
tween individual spins. This suggests that the QCG de- 
localization transition is actually in the same universality 
class as the MBL transition. The above scenario con¬ 
trasts an alternative picture of the QCG-delocalization 
transition as being driven chunks of perfect metal inter¬ 
acting through QCG regions via stretched-exponential- 
in-distance couplings, which would predict an absence of 
snbdiffusion in the delocalized phase |l9| . The above con¬ 
siderations instead suggest that the “metallic” clusters 
in this picture are better thought of as “critical” clus¬ 
ters, and interact with highly renormalized exponential- 
in-distance interactions. 

We note that similar arguments can be used to sug¬ 
gest the stability of QCGs with stretched exponential 
interactions. MBL phases with (short-ranged) power-law 
interactions, ^(a:) ~ and MBL in higher dimensions 
for strong disorder. Loosely, though the two-spin inter¬ 
action. V(L) between spins of scale L, is much larger 
than the many-body level spacing 5{L) ~ 2~^ (where d 
is the number of spatial dimensions) for large L. in or¬ 
der to access the small many-body level spacings, ~ 
powers of the two-spin interactions must be used. For 
strong disorder, each of two body transition is energy- 
mismatched sucL that the many-body matrix element, 
r(L) corresponding to transitions with level spacing 5{L) 

also scales as r(L) ~ (^) allowing for an MBL phase 
for W ^ V. Similar arguments to those above suggest 
that the MBL transition with extended interactions will 
be of the same imiversality class as that with strictly 
short-ranged interactions - due to the renormalization of 
the interactions between resonant clusters. We leave a 
more detailed study of these issues for futine work. 


C. Validity of quantum percolation model 

Before concluding, we pause to comment on the as¬ 
sumptions and approximations that have gone into con¬ 
structing the effective quantum percolation model of the 
delocalization transition. Our key assumption is that 
the critical resonant cluster has a self-similar structure 
in energy and space, which allows us to hierarchically 
construct the percolating resonance that drives delocal¬ 
ization. This hierarchical construction cleai'ly does not 
describe the most general collective resonances, and in 
particailar ignores generic unstructured resonant clus¬ 
ters that do not decompose into a hierarchical tree of 
few-body resonances. However, the assumption of self¬ 
similarity is naturally expected for a continuous phase 


transition characterized by a diverging length scale, 
based on knowledge of more conventional disordered crit- 
icahty. We therefore expect our approach to accurately 
capture the universal aspects of such a resonant structure 
close to the delocalization transition. 

Our approach also seeks to construct a delocalizing res¬ 
onance by merging metallic' clusters. Such an approach 
is known to be problematic for non-interacting Ander¬ 
son insulators in which the conductivity of blocks does 
not add in a simple fashion due to coherent interference 
effects [35]. However, in the present context, since the 
system must be thermal at the MBL transition, we ex¬ 
pect strong dephasing to wash out any effects of quantum 
interference (e.g. the time-scale for dephasing is more 
rapid than that of transport, and does not diverge even 
in the MBL phase). 

VII. DISCUSSION 

To summarize, we have developed an effective model 
for the many-body delocalization transition of a one¬ 
dimensional MBL phase. This model can be efficiently 
simulated for very large system sizes (with computational 
cost scaling polynomially in system size), and we argue 
that it gives an accurate description of its long-distance, 
low-energy, universal critical properties. 

This transition represents a novel form of non¬ 
equilibrium excited-state criticahty, distinct from clas¬ 
sical, thermal phase transitions and zero-temperature 
quantum-critic'ality. We find that this excited state tran¬ 
sition is continuous (second order), characterized by a 
single diverging length scale. We describe the set of uni¬ 
versal critical exponents that characterize this melting 
transition, which are connected by scahng relations. In 
particular, we find an muisual scaling structure of the 
critical slowing down of dynamics and thermal transport 
characterized by a continuously evolving dynamical ex¬ 
ponent, z, which diverges in a imiversal fashion upon 
approaching the quantum critical point. This unusual, 
universal divergence of a critical exponent, z, stems from 
a sparse structure to the transport path tlnough the crit¬ 
ical percolating cluster. It is natural to expect that this 
sparsity corresponds to a (multi)-fractal structure of the 
critical cduster. While our approach makes many simpli¬ 
fying approximations, and is clearly not exact, we believe 
the general scaling structiue it predicts (Fig. to be cor¬ 
rect, and moreover expect to obtain reasonable (though 
possibly approximate) value for the correlation length ex¬ 
ponent, u. In particular, it is encouraging that the scaling 
data of Fig. I^is closely compatible with the renormaliza¬ 
tion group approach of Ref. [15] |33|, despite the marked 
differences between these approaches. Testing these re¬ 
sults using standard mmieric'al techniques may be chal¬ 
lenging. In particular, a direct frontal assault via exact 
diagonalization is unlikely to give accurate scaling results 
in systems of only a few tens of sites. 

We anticipate that these methods can be extended 
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to treat a variety of other out-of-equilibrium dynamical 
phase transitions in disordered quantum systems. For 
instance, our analysis in terms of quantum percolation 
of resonant clusters can be straightforwardly adapted to 
treat higher dimensional delocalization transitions, which 
are so far inaccessible by other methods. One also ex¬ 
pects a self-similar fractal structure of the critical res¬ 
onant cluster, which may have implications for certain 
scaling properties and many body level statistics at the 
delocalization transition. A key expected difference in 
higher dimension is that the near-critical transport will 
not be dominated by rare insulating gaps, which can be 
simply circumvented in higher dimensions. This sug¬ 
gests that the delocalized thermal liquid will be diffu¬ 
sive with dynamical exponent z = 2. Then repeating the 
argnments leading to the scaling relation C, = i-' above, 
one would expect the diffusion constant to vanish as 
D ~ upon approaching the higher-dimensional 

MBL transition. 

Another unanswered question regarding one¬ 
dimensional systems concerns the fate of subdiffusion 


deep in the delocalized phase at weak disorder, far from 
the MBL transition. One possible scenario is that sub¬ 
diffusion persists to arbitrarily weak disorder, such that 
diffusion is only asymptotically recovered in the clean 
limit. This would imply that any amount of disorder in 
one-dimensional systems leads to subdiffusive dynamics. 
The alternative would a second finite-disorder-strength 
excited state phase transition separating diffusive and 
subdiffusive regimes. We leave a detailed investigation 
of such issues for future work. 
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Appendix A: Details of Effective Model of the 
Delocalization Transition 


1. Resonance Merging Procedure 


In the first stage of resonance formation, we check for 
pairwise resonances between two spin pairs i and j for 
which r(a:ij) > SEij, where r(a;) ss is the typi¬ 

cal effective coupling through the intervening strong dis¬ 
order phase given by (j^. Here, Xij is the distance sepa¬ 
rating the spin pairs i,j and 5Eij is their energy differ¬ 
ence. As we iterate this procednre we generate larger 
clusters and check whether excitations can resonantly 
tumiel between any two clusters. Specifically, at any 
stage in the iterative process, we have a collection of N 
resonant chains labeled hy j = 1... N, each comprised 
of [rrij + 1) bonds. Each resonant cluster has ss 2’”-> en¬ 
ergy levels, spread over bandwidth Aj equal to the sum 
of couplings ,7 for the resonant links inside the cliLSter, 
and hence have typical level spacing 5a ^ 2 ^' 
ergy levels within a given resonant chain are all strongly 
a<lmixed and exhibit level repulsion. 

If the residual coupling r(a;ij) between two clusters i, j 
can resonantly drive transitions among the levels, then 
we merge these two resonances uito a new, larger one. 
The precise rule for whether or not to merge two res¬ 
onances depends on the relative size of the bandwidths 
and level spacings for each. Consider two resonant clus¬ 
ters, labeled j = 1.2. Without loss of generality, assume 
that cluster 1 has larger bandwidth Ai > A 2 . There 
are two distinct cases to consider. If the level spacing 
of the first cluster exceeds the bandwidth of the second, 

> A 2 (Fig. Cl i), then the minimum energy difference 


(a) 



Case 1; (5i > A 2 


(6) Case 2; < A 2 


^■ 1 = 4 .— 

Ai - + .-I — =i> 

I — t.IiL- 


Cliister [1] CliLSter [2] 



[ 1 ] ® [ 2 ] 


FIG. Cl. Two cases for merging resonances. The energy 
cost for changing the state of one cluster (blue circle to open 
circle connected by black arrow) can be partially compensated 
by changing the state of another cluster. The reduced energy 
penalty, SE 12 , for such “flip-flop" transitions is either Ji — A 2 
or depending on whether 5i ^ A 2 , as shown in panels 
(a) and (b) respectively. 


associated with changing the states of both resonances 
is; (5 £* 1,2 = <5i — A 2 . On the other hand, if Si < A 2 
(Fig.lClb), then the minimum energy cost for transitions 
is of order (5£’i,2 = 

If r(xi 2 ) ‘C then the residual coupling between 

the two resonant clusters merely weakly perturbs the en¬ 
ergy levels of each. In this case, the combined energy 
levels of both clusters are well-described by the tensor 
product of the independent levels associated with each 
cluster indicating that energy excitations are well local¬ 
ized independently in each cluster. On the other hand, 
if r(xi 2 ) ^ 5Ei^2, then the energy levels of both res¬ 
onances are strongly admixed such that excitations are 
spread throughout both clusters, and energy can reso¬ 
nantly shuttle between them. In om numerical proce- 
diue, we take 6E as a sharp cutoff, merging resonances 
coupled by r(xjj) > SEij. 


2. Renormalization of inter-cluster couplings 

In this Appendix we explain the conceptual necessity 
of renormalizing the inter-cluster couplings in the effec¬ 
tive resonance percolation model for the MBL transition. 
Perhaps surprisingly, we find that even though the renor¬ 
malization steps are conceptually necessary to obtain a 
well-defined transition in asymptotically large systems, in 
practice, they may be omitted without changing the uni¬ 
versal scaling properties obtained from finite-size scahng 
in moderately large systems. 

If one omits the renormalization of inter-cluster cou¬ 
plings, and simply takes bare couplings; F^- = , 

where Xij is the distance between clusters i and j, there 
is. in principle, a rare event chiven breakdown of the 
strong-disorder phase which rounds the transition in the 
effective resonance percolation model [fl.*!! 124) . Namely, 
consider the case where one finds a resonant cluster of 
N, spins, such that the bandwidth of this cluster An. is 
larger than the maximal energy mismatch between two 
spins (« W). Then, with simple geometrically dictated 
couplings that depend only on the inter-spin distance, 
this Nq size cluster will deterministically absorb its neigh¬ 
bors, and continue to grow in this fashion until it encom¬ 
passes all spins. Such an “avalanche” style breakdown 
occurs even for arbitrarily large disorder, as there is al¬ 
ways a non-zero probability of finding such an N, spin 
“seed” cluster. 

Specifically, for very strong disorder, ^ 1 

spins are required to form a cluster whose bandwidth. 
Fat. ^ N^V, exceeds the energy mismatch W. The prob¬ 
ability of finding such an iV,(l'F) spin resonant cluster is 
Favaianche ~ , indicating that an avalanche will oc¬ 
cur for systems larger than Tavaianchi- ■ While this 

system size is extremely large for strong disorder W ^V, 
it nevertheless indicates that the non-renormalized cou¬ 
plings lead to the unphysical prediction that the MBL 
phase will always be unstable to such an avalanche-style 
breakdown, and shows that the inter-cluster renormal- 
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FIG. C2. Finite-size scaling for the disorder averaged 
localization length fjoc (See Fig.^ using a simplified model 
that ignores the renormalization oithc couplings. Results are 
averaged over ~ 10'* disorder realizations. 


ization steps are necessary to obtain an asymptotically 
well-defined phase transition. This avalanche instability 
is eliminated by proper renormalization of inter-cluster 
couplings. Again working in the limit of strong disorder, 
we observe that, as the “seed” cluster is formed, its cou¬ 
plings to neighboring spins will be typically renormalized 
downward by a factor of ^ for each spin in the “seed” 
cluster. Hence bai'ring an extensive set of accidental reso¬ 
nances. which has vanishing probability for large systems, 
the avalanche will be suppressed for ^ <C 1. 

On the other hand, simply omitting the renormahza- 
tion steps and proceeding with the bare couplings, we 
obtain essentially identical finite-size scaling results com¬ 
pared to the case with renormalized couplings (Fig. C2 i. 
In particular we obtain finite size scaling curves for 
that cross nicely at an apparent critical disorder strength 
^ R3 17 with V = 0.3. While the precise (non univer¬ 
sal) value of \Vc differs in the absence of renormalization 
steps, the universal scaling structure such as the crossing 

S3 2 • 10“^ and correlation length exponent 

w, 

u w 3.5 i 0.4 agree witliin error bars with those results 
presented in the main text. 

Thus, the theoretically predicted avalanche phenom¬ 
ena appears to be completely absent even for relatively 
large (L ss IQ^) system sizes. In fact, at the MBL tran¬ 
sition point ^ ss 17. from the arguments above, the 
avalanche instability would only be visible for system 
sizes of L > spins, wliich is completely unaccessible 
by either numerical simulations or, for that matter, cold 
atoms experiments. Hence, while the renormalized cou¬ 
plings are in principle crucial for obtaining a well-defined 
transition, in practice one may omit the computationally 
expensive renormalization step and still observe miiver- 
sal critical finite-size scaling behavior in large systems. 


value ^ 


This situation is somewhat analogous to accessing criti¬ 
cal properties of a conventional quantum critical point at 
very small, but finite temperature T ss , which 

in cuts off the flow towards the imiversal scaling regime, 
but only above a very large length-scale, helow which 
is it possible to accurately identify the universal scaling 
exponents. 


3. Construction of transport path 

Since we have seen that the potential avalanche in¬ 
stability is absent for computationally accessible system 
sizes at the critical point (W = we omit these 

computationally expensive steps in computing transport 
properties. However, in this caise, some care is required to 
properly construct the transport path through the reso¬ 
nant cluster, and to avoid a related avalanche-style break¬ 
down of the transport calculation for W < W^. Namely, 
despite the absorption of all bonds into the resonant clus¬ 
ter for W < Wc, not all bonds contribute equally to 
transport, and some care is required to properly iden¬ 
tify the transport path thi'ough the resonant cluster. An 
illustrative example is shown in Fig. |C3| 

At some intermediate step, i — 1, in the resonance¬ 
merging process small resonantly linked clusters of spins 
(blue rectangles) of various sizes have been identified. At 
step i, smaller clusters merge to form two large reso¬ 
nantly linked clusters, Ci,C 2 . In the next step, i -I- 1, 
after accounting for the reduced level spacing of the 
newly formed clusters, one may find that Ci and C 2 can 
resonantly exchange energy via the black link, and are 
merged into a single larger cluster Ci U C 2 - It frequently 
happens that, at step i + 1, the level spacing of Ci (or 



FIG. C3. Transport path through resonant cluster 
(Ik < Wc)- Schematic depiction of sequence of resonance- 
merging steps (i. i -h 1, i + 2) that produces a gap in the trans¬ 
port path through the resonant cluster. See text for detailed 
explanation. Blue rectangles are resonantly linked clusters of 
spins. Black and gray lines indicate resonant links between 
previously formed clusters. Dashed crossed out red line indi¬ 
cates a spurious connection that should not be added to the 
transport path of the cluster. 













12 


C 2 ) alone is too large to facilitate the resonant exchange 
of energy between smaller clusters, e.g. C 3 , between Ci 
(or C 2 ), but that the combined level spacing of Ci U C 2 
is small enough to resonate with C 3 . Then, in the sub¬ 
sequent step, 1-1-2, Cl U C 2 will absorb the intervening 
spins. However, it is crucial to note that, since C 3 can¬ 
not (in this example) resonantly exchange energy with 
Cl or C 2 alone, but rather only by collectively exciting 
both Cl and C 2 . Hence, resonantly transferring energy 
between C 3 and Ci U C 2 takes at least as long as the 


time. t{Li 2 ), to tunnel between Ci C 2 (long black 
link). A computationally inexpensive way to account for 
this feature, is to draw only a single connection between 
the closest pair of bonds on C 3 and Ci U C 2 (gray lines 
ill Fig. |C3[ l. Drawing additional connections, such as 
the crossed-out dashed red line in Fig. |C3| would spu¬ 
riously allow energy to hop from Ci —)• C 3 —^ C 2 in 
time t(Li 3 ) + r(L 32 ) ‘C r(Li 2 ), whereas in this example 
we have seen that transferring an excitation from C 3 to 
Cl U C 2 must take at least t(Li 2 ). 


